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Constitutive equations are proposed for dense materials, based on the identification of two types 
of free-volume activated rearrangements associated to shear and compaction. Two situations are 
studied: the case of an amorphous solid in a stress-strain test, and the case of a lubricant in tribology 
test. Varying parameters, strain softening, shear thinning, and stick-slip motion can be observed. 
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From food to beauty products, from our bone joints 
to the gouge of tectonic faults, the dynamical properties 
of dense materials determine important and ubiquitous 
phenomena that govern our life. The idea emerged re- 
cently, that some sort of universality might be at work 
in structural systems, ranging from glasses to granular 
materials. Q In this letter, I hypothesize that univer- 
sality results from the crucial role played by structural 
rearrangements. The approach I propose relies on the 
so-called shear transformation zone (STZ) theory, intro- 
duced to account for elasto-plastic transitions in amor- 
phous solids. Q STZ theory provides a scheme for re- 
arrangement kinetics in a dense material: Macroscopic 
deformation is seen as result of local free-volume acti- 
vated rearrangements at a mesoscopic scale. Local states 
are introduced, related to the local orientations of con- 
tact network. Mean-field equations of motion for those 
local states lead to macroscopic constitutive equations. 

STZ theory is originally designed for solids, (|] but 
it has been shown recently that it could be adapted to 
the case of granular materials and successfully account 
for important features of granular flows. [^) It is shown 
here that complementing STZ theory with free-volume 
kinetics permits to account for important features of 
both amorphous solids (strain softening) and dense liq- 
uids (stick-slip) within the same framework. In previous 
works on STZ theory free-volume appears as a parame- 
ter, although it is obviously a state variable that varies as 
the system dilates or contracts. Free-volume kinetics is 
inspired by recent results on granular materials submit- 
ted to vertical tapping: local rearrangements associated 
to a buckling instability of chain forces have been shown 
to account for logarithmic relaxation of free-volume. ||] 
The micro-structural picture proposed for dense materi- 
als involves two types of rearrangements, associated re- 
spectively to shear motion and to free- volume relaxation. 
The interplay between those two mesoscopic processes 
and associated timescales is shown to account for impor- 
tant properties of sheared materials. 

To emphasize the analogy between solids and dense 
liquids, two important experimental tests are considered: 
the case of a stress-strain test performed on a plastic 



material ||] and the case of an (overdamped) tribol- 
ogy test where a dense liquid or a granular material is 
sheared. ||-|| In the former case, the resulting constitu- 
tive equations are shown to account for strain softening. 
In the latter, they account for a transition between steady 
sliding and stick-slip motion. The picture is completed 
with the case of constant stress experiments (creep tests) . 

Let me now introduce a basic mechanism for rearrange- 
ment kinetics, along the lines of STZ theory. ^ The 
present approach remains at a mean-field level. The over- 
all shear deformation is denoted e and decomposed into 
an elastic and a plastic part: e = e el +e pl . The elastic de- 
formation is proportional to the shear stress a = fie el ; the 
plastic deformation results from local, irreversible trans- 
formations of the contact network. A local rearrange- 
ment involves several molecules at a mesoscopic level. A 
shear transformation zone is defined as a locus within the 
material where an elementary shear is made possible by 
the local conformation of neighboring molecules. An es- 
sential remark that lies at the root of STZ theory is that 
once a rearrangement has occurred, some contacts break, 
some others are formed, and the molecules involved can- 
not shear further in the same "direction" , although they 
might shear backward. This leads to identifying pairs 
of types of arrangements, that are transformed into one 
another by a local, elementary shear. To simplify the the- 
ory, one pair of orientations is considered, corresponding 
to the principal axes of the stress tensor. An elementary 
transformation can be sketched as follows, 




although an actual rearrangement involves in general 
more than four grains. Populations (number densities) of 
arrangements of each type are denoted n± , and R± de- 
note the rates at which a ± — > =p transformation occurs. 
The rates R± depend on force and free-volume fluctua- 
tions and will be defined further. The macroscopic plastic 
shear is evaluated from the balance between both types 
of elementary motion: 

e pl = Aq {R+n + - i?_n_) , (1) 
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where some constant Aq has been introduced. Equa- 
tion ([!]) must be complemented by equations of motion 
for the populations n±. They are written of the form, 0] 



n ± = R T i2 T — R±n± + a i pl (A c — A a n±) 



(2) 



A c and A a are constants. The first two terms on the rhs 
account for the internal reconformations of STZ's, while 
the last term introduces a coupling with the mean flow: 
At the macroscopic scale, the flow constantly stirs the 
molecules, thus creating and destroying local configura- 
tions. The rate at which the macroscopic flow induces 
new configurations is estimated as the ratio of the over- 
all work, a e pl , over by some (constant) typical normal 
force (the force required to modify a contact). 

Transformation rates R± depend on volume and force 
fluctuations. A local reconformation occurs if some vol- 
ume is available around a transformation zone, and if 
the local force network is appropriately oriented. In this 
work, I assume that the probabilities associated with 
those two types of fluctuations factorize. This differs 
from the prescription given in || . An exponential distri- 
bution of volume fluctuations is assumed, and vo denotes 
the typical free- volume per grain needed for a rearrange- 
ment. The probability that sufficient free- volume is local- 
ized around a given STZ is proportional to exp (—vq/v/), 
where Vf is the average free- volume. The distribution of 
forces is also expected to be exponential. Shear stress 
a introduces a bias of the force network. The probabil- 
ity that some strong local bias of the force network is 
realized and promotes shearing in the ± — > =p direction 
is estimated to be proportional to exp(±<r//Z), where j2 
measures some typical force of the stress network. In the 
following, = 1 is imposed, which fixes the unit offerees. 
The rates are thus written, 

R± = Rq exp (-v /vf) exp (±a) , 

with Rq, the update frequency of microscopic processes. 
This choice of an expression for the rates R± will be dis- 
cussed further, but other choices are possible that lead 
to similar results. M 

At this point, free-volume appears as a parameter. 
However, this picture is now complemented with free- 
volume kinetics. Free-volume relaxation results from 
tiny reconformations of molecules into better local pack- 
ings, pj A second type of rearrangement is thus intro- 
duced, which is essentially isotropic — hence does not 
require any particular oricntational zone to occur — but 
which is otherwise treated in a similar and consistent way. 
The typical volume fluctuation needed for a local collapse 
is denoted v\ ; free- volume relaxation occurs at a rate pro- 
portional to exp(— vi/vf). Note that v\ is expected to be 
smaller than vq as an improvement of the local packing 
should require reconformations of a smaller scale than 
shear motion. In the absence of shear (a = 0), and at 
temperature, the equation of motion for free-volume is 



expected to reduce to, if = —E±exp(—Vi/vf), which 
accounts for a slow (logarithmic) relaxation of v/. Q] 
When the material shears, some dilatancy is expected. 
This is due to the fact that macroscopic shear motion 
brings about new random configurations of molecules. 
In new arrangements, the local free-volume is not opti- 
mized. This feature must be incorporate in free-volume 
dynamics via a creation term. To estimate this creation 
term, the following argument is proposed: the average 
dilatancy due to shear deformation is viewed as a trans- 
fer of a fraction of the energy aSc pl into PSvf by some 
sort of levering effect within the contact network. Pres- 
sure dependence is incorporated in constants. When the 
system deforms by 5e pl , the average variation Svf of the 
free- volume is thus expected to be proportional to aSe pl . 
Equation of motion for Vf results from the interplay be- 
tween local relaxation and shear-induced dilatancy. It is 
written: 



Vf 



-E 1 exp [—vi/vf] + A v cre 



pi 



(3) 



In the following, constant A v = 1 is imposed, which fixes 
the unit of Vf . 

Before proceeding further, equations (|l]) and (||) are 
written in a more appealing form. Following 0J, vari- 
ables 



A = 



n_ — n + 



and, A = 



n_L + n_ 



are introduced, along with the rescaled parameters — 
2A C /A a , e = AoA c /Aa, 7 = A) Ac, and E = 2e a R . 
It comes, 



e pl = E exp [-v /vf] (A sinh(cr) - A cosh(cr)) 
■pi 

A = — (1 - 7<r A) 
■pi 

A = 7cr— (1-A) 



(4) 
(5) 

(6) 



Variables A and A represent respectively the total nor- 
malized density of STZ's and the bias between the pop- 
ulations n±. In the absence of shear motion, the dy- 
namics of A and A freeze, while Vf undergoes some "au- 
tonomous" relaxation: those two sets of internal variables 
are associated to very different types of memory. The ini- 
tial values given to those variables account for the way 
the system has been prepared. 

Let me comment on initial values used for free- volume. 
In the absence of shear, a = 0, the system still undergoes 
free-volume relaxation. If the picture of a dense mate- 
rial was complete, a reasonable initial state should cor- 
respond to values of Vf for which the factor exp(— vi/vf) 
(hence exp(— v /vf)) is vanishingly small: the system is 
an aged glass. However, for a finite temperature, free- 
volume is expected to equilibrate at some non-vanishing 
value. Temperature dependence is not included in the 
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current approach. In this letter, to circumvent this ques- 
tion, the initial condition for vt is a fraction of vq, which 
is small, but not as small as would be expected from an 
aged glassy material. In fact, for very small initial values 
of the free-volume, a large stress is necessary to initiate 
deformation. This produces behavior analogous to the 
so-called memory peak. Those aspect are not considered 
here. 

I first describe briefly the resulting behavior in a con- 
stant stress experiment (creep test), before discussing the 
case of constant strain-rate and tribology experiment. 

Stress a is fixed. To simplify the analysis, I separate 
equation from equation (||). Equations (^-Bl) admit 
two types of solutions: multiple jammed states 
provided A = tanh(er)A, and a flowing regime 
with A = 1 and A = I/7CT, whence 



■pi 



0), 
^0), 



■pi _ 



E exp[— vq / v f](sinh(a) — cosh(a) / (ja)) . 



The jammed states are stable for small values of the ap- 
plied stress, and the two types of solutions exchange sta- 
bility when the yield stress a y is reached, with a y so- 
lution of crtanh(cr) = I/7. (2) The rhs of equation (||) 
is always negative if a < a c , where a c is the solution 
of E (a sinh(c) — cosh(cr)/7) = E\. In this case, the 
system remains "solid", free- volume relaxes, hence e p : 
the material creeps. Above ct c , equation (|^) admits one 
unstable fixed point: depending on its initial value, Vf 
either decays or diverges. In this latter case, the weights 
exp(— Vi/vf) and exp(— vo/vf) saturate to 1, and e pl goes 
to a constant. The latest stages of the dynamics, when 
v f gets large, are beyond the scope of the current work. 
Various creep curves are shown figure p\- (left) which com- 
pare with experiments presented in 
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FIG. 1. Left: Creep test, deformation as a function of time 
at fixed a. Parameters Eo = E\ = 0.1, vo = 0.5, «i = 0.4, 
7 = 10, eo = 1, and initial conditions A = 0, A = 0.01, and 
v f = 0.25. From bottom to top, a = 2, 3, 3.7, 3.9, 4.1, 4.2, 
4.5. Right: Friction force (oc a) as a function of the velocity 
of the upper plate (oc e pl ) during stick-slip motion as shown 
figure I (right). 



In a stress-strain test, 

c el ' 



e is fixed, 
or a = n(e 



The shear stress is 
determined by a = ^te e( , or a — /i(e — e pl ). Let me show 
that the same relation appears for overdamped stick-slip. 
A layer of liquid is sheared between two plates. The gap 



between the plates is denoted a, and the surface of con- 
tact S. A pulling force is exerted on the upper plate, 
with a spring of stiffness k, pulled at velocity V. If the 
inertia of the upper plate is neglected: F — k(Vt — x). 
The position x of the upper plate is related to the defor- 
mation of the liquid by: x = 2ae pl , and the shear stress is 
a = F/S, whence, a = {k/S){V - 2ae p/ ), which is of the 
form & = /.i(e — e pl ). Identical equations govern both sys- 
tems, although, remarkably, a solid is usually considered 
in the former case, and a liquid in the latter. 

In the following, A = 1 is taken as initial condi- 
tion. This is consistent with the fact that stress-strain 
curves (or stiction peaks) present well-defined features 
after preparation of the sample. Since A saturate to 1 
after some shearing, it cannot govern some repeatable 
test, but accounts for the evolution of the material dur- 
ing preparation. 
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FIG. 2. Behavior of sheared materials for fixed strain rate 
e. Parameters of the theory are, Eo = 0.1, E\ = 0.1, vo = 0.5, 
v\ — 0.4, 7 = 10, eo = 1- Initial conditions, are A = 0, A = 1, 
Vf = 0.25. Left: strain softening, for e = 0.08 (solid line), 
0.1 (dashed line), 0.15 (dot-dashed line), 0.2 (dotted line). 
(Those values are defined up to a time constant, and do not 
need to compare quantitatively with actual values.) Diver- 
gence of Vf is seen for the three largest values of e. Right: 
stick-slip motion obtained for e = 0.01. 

Example of behavior resulting from equations (|[^) is 
shown figure ^. For the large values of e, the system fol- 
lows a stress-strain curve which present strain-softening. 
(Note that this is for the chosen of parameters, and that 
strain hardening can be seen for other values.) In the lan- 
guage of tribology, a stiction peak leads to steady slid- 
ing. Softening can be understood as follows: During a 
first stage of the dynamics, plastic deformation e pl re- 
mains zero, or very small, either because, a is below the 
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yield stress, or below the free- volume remains small. The 
shear force increasing, transient creeping motion induces 
dilatancy, thus triggering plastic deformation, further in- 
crease of Vf, and a starts to decrease with the increas- 
ing activation factor exp(— vo/vf). For the largest values 
of e (here, e = 0.1,0.15,0.2), the free-volume diverges, 
leading the system to melting or break-up. For inter- 
mediate values (here, e = 0.08) the system converges to 
steady sliding with finite free-volume. For small e, stick- 
slip motion is observed, sliding peaks being associated 
to sudden rises of the free-volume. The friction force is 
multi- valued, and a as a function of e pl is shown figure [j] 
(left). The cycle is followed in the clockwise direction and 
is remarkably similar to cycles reported in experiments 
on lubricants. 

Let me now comment further on the expression cho- 
sen for the transformation rates R±. Other prescriptions 
can be used, [§] but introduce some correlations between 
free- volume and force fluctuations. It has been checked 
that the results presented here do not depend on other 
possible choices of R±, so long as free- volume activation 
is incorporated. The factor exp [—vo/vf] is essential, and 
accounts for the fact that at low free- volume, activated 
processes slow down dramatically. The preceding expres- 
sion is chosen because it allows to factorise exp [—vo/vf] 
in equations (|l|-^|). It simplifies somehow the resulting 
constitutive equations, but more importantly, the free- 
volume enters the kinetics of shear motion only as much 
as it determines its time-scale. The coupling between 
free- volume relaxation and shear deformation is realized 
through this mechanism only, which is a rather minimal 
and subtle effect. It is profound and particularly illumi- 
nating that such a subtle coupling via time-scales suffices 
to obtain so many non-trivial features. 

It is noteworthy that early approach to stick-slip relied 
on the introduction of rate-and-state laws, in the spirit of 
friction laws used to study earthquakes dynamics. fic| l In 
the current approach, free-volume is identified as a loga- 
rithmically relaxing state variable responsible for strain- 
softening and stick-slip motion. The importance given to 
free- volume provides a direct test of the theory since some 
appropriate control of the free-volume kinetics could, in 
principle, prevent stick-slip motion. Although difficult to 
realize in most experimental set-ups, this might turn out 
to have an important practical impact. 

The current approach has remained at a qualitative 
and general level. This theory has now to pass the test 
of more stringent comparison with experimental data. 
This might require minor adaptations, e.g. to the case 
of granular materials, ||,^| but the mechanism proposed 
is of great simplicity and should adapt to various types 
of structural systems. A complete picture of dynami- 
cal properties of dense materials will require, however, 



to incorporate spatial extension of a sample of material 
and account for localization of the deformation into shear 
bands. Recent developments in STZ theory have shown 
a possible mechanism JTT] | and should benefit from the 
incorporation of free- volume dynamics. 
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